Signatures of electron correlations in the transport properties of quantum dots. 
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The transition matrix elements between the correlated N 
and N+l electron states of a quantum dot are calculated by 
numerical diagonalization. They are the central ingredient 
for the linear and non-linear transport properties which we 
compute using a rate equation. The experimentally observed 
variations in the heights of the linear conductance peaks can 
be explained. The knowledge of the matrix elements as well 
as the stationary populations of the states allows to assign the 
features observed in the non-linear transport spectroscopy to 
certain transition and contains valuable information about 
the correlated electron states. 
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By using modern nanostructure fabrication technology 
a few electrons can be confined to very small regions 
in space [Q. In these so-called quantum dots or arti- 
ficial atoms the Coulomb interaction between the elec- 
trons is very important for understanding their quantum 
mechanical properties. Weak coupling to external reser- 
voirs via tunnel barriers allows to observe single electron 
transport effects like the Coulomb blockade oscillations in 
the linear conductance at millikelvin temperatures ■ 
In non-linear transport, features are observed which are 
closely related to the excitation spectrum of the interact- 
ing electrons S. 

Transport involves transitions between the many-body 
eigenstates of the confined electrons. They are approxi- 
mated as products of one particle states within the charg- 
ing model where the Coulomb interaction is modelled 
phenomenologically by the capacity of the quantum dot 
||. This is not sufficient to explain e.g. the experimen- 
tally observed negative differential conductances M . Es- 
pecially at low electron densities, correlations between 
the electrons are crucial. This was explicitly shown for 
a quasi one-dimensional (ID) box || where the correla- 
tion leads to TV pronounced peaks in the charge density 
distribution if the mean electron distance r s exceeds the 
effective Bohr radius with = [m e / m)eas (e rela- 
tive dielectric constant, m effective mass). In the present 
paper we use the same model to investigate the influence 
of the spatial properties of the correlated many-electron 
states on the linear and non-linear transport properties 
of an artificial atom. An investigation in a similar spirit 



was recently performed for a harmonic confining poten- 
tial in 2D However transport properties were not 
calculated and only the transition between N = 2 and 
3 electron states was studied in this work. We demon- 
strate here that the current-voltage characteristics ob- 
tained by solving a stationary master equation exhibits 
very specific signatures of the electron correlations. They 
influence the transition matrix elements and also station- 
ary occupation probabilities of the states. Some of the 
" lines" observed in non-linear transport spectra are even 
enhanced. Such a result cannot be obtained by consider- 
ing only the transitions rates. 

We will show that apart from the restrictions due to 
spin selection rules discussed earlier JuJjnJ transitions 
are suppressed or enhanced when taking into account the 
spatial properties of the wave-functions. The heights of 
the peaks in the linear conductance become non-equal 
even without taking into account the energy dependence 
of tunneling matrix elements. In non-linear transport 
the excited levels of the quantum dot which can be ob- 
served are closely related to the most prominent allowed 
transitions and to the highest stationary occupation prob- 
abilities. 

As a model for the quantum dot we consider a quasi- 
1D square well §,0 of length L = 9.5a* and N < 4. 
The corresponding mean electron density is close to the 
one in experiments on GaAs-AlGaAs heterostructures 
where the mean distance between the electrons is about 
3a* (wlOnm) |gg. 

Including the spin degree of freedom a, the dot Hamil- 
tonian reads 

H D = X]( £ » ~ e( t>) C n,o C n, a 

+ Kl4K37l2ni C n4CT C i 3 <T' C ™20-' C «lO- ■ W 



The electrostatic potential <j> depends on the gate and 
transport (bias) voltages which are applied to the sys- 
tem and shifts the energies of the one-electron lev- 
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is the matrix element of the interaction 
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/ey/(x — x 1 ) 2 + A 2 . The cutoff at short dis- 
tances simulates a small transversal spread of the wave 
functions (A<I) and provides finite K l4n3 „ 2ni . 

In the limit of high tunnel barriers, transport is de- 
termined by the eigenstates of the isolated dot. From 
the latter we calculate the transition probabilities be- 
tween N and (N ± 1) electron states. Each energy 
eigenstate is simultaneously eigenstate of the total spin 



S* 2 = (<5"i H h <tn) 2 which implies a (25+l)-fold de- 
generacy with respect to S z in the absence of a magnetic 
field. 

We include n= 1, . . . , M one-electron states <£n(a;)x<T 
when diagonalizing i/rj. Here, tp n is a spatial one electron 
function and \„ a spinor with er =J,, t • The Hamiltonian 
matrix in the basis of Slater determinants is of the rank 
of the binomial number r — C^ M (r < 1.5 • 10 4 , for 
M = 10, . . . , 13 ). The Lanczos method was used when 
r > 3 • 10 3 . 

By its algorithm |l3| the Lanczos diagonalization pro- 
vides only one eigenvector | 'I' 5 ) Lane for each energy eigen- 
value. The calculation of transition rates is consider- 
ably facilitated when using eigenstates of S z . Usually 
is a linear combination of all of the 2S+1 vec- 
tors in the subspace of Zeeman levels. In order to recover 
the eigenvectors of S z we apply projectors Ps z IvE^Lanc °c 
\E ,S ' -Sl } corresponding to a specific S z . After normaliza- 
tion IvJ/ 5 -^) = J2v=i bf' Sz W) can be expanded into the 
Slater determinants \v) = cL ff • • • c\ lN<7N |0) of the non- 
interacting electrons. The coefficients b„' * are obtained 
after diagonalization, projection and renormalization. 
By construction we have bf = if o\ + • • • + o ^ ^ S z . 
The method would fail in the unlikely case that )Lanc 
is accidentally perpendicular (within the numerical accu- 
racy) to one of the |\E ,S '' S *). The procedure can also be 
applied to higher dimensional models. 

To study transport properties, we use the tunnel- 
ing Hamiltonian and the usual rate equation approach 
|uj[nj. Then the dc-current through the quantum dot 



I L/R = (-/+)e P Kj R ( N i 



Ni) (2) 



i,3 U¥=i) 



corresponds to the rate of electron passages through the 
left or the right barrier. I is computed from the station- 
ary occupation probabilities Pj which are solutions of the 
equation ^ C^i jPj — Tj jPj) = 0. The transition rates 

between all of the many-electron states indexed by j are 

r, 



3,1 



r L . 



r R 



Time dependent perturbation theory yields 



^L/R 



t L/R 



,L/R 



i\ c n,& 



f L/R (E) 



[1-/ L/R (-B)], (3) 



for transition probabilities between the eigenstates 
and | tyj ) of Hu with Nj — Ni + 1 in lowest order in the 
tunneling. The indices i and j contain in particular the 
electron number N, the total spin S and S z . t L / R are 
the tunneling rates through the left/right barrier. The 
electron has to provide the energy difference E = Ej — Ei 



when entering or leaving the dot. The Fermi-Dirac dis- 
tribution functions f L ^ R (E) describe the left/right reser- 
voirs with chemical potentials /i L / R . 

The energy spectrum of the N electrons for the den- 
sities studied consists of multiplets [p||l^]. The energy 
differences fi between the latter are considerably larger 
than the intra- multiplet energy differences A. They are 
important only for large transport voltages. For a GaAs- 
AlGaAs heterostructure of length L — 9.5a# with N = 4, 
£1 fss Q.2meV, and A ks 62fieV. For small transport 
voltages we can restrict ourselves to transitions between 
states within the lowest multiplets (table |). The total 
number of states within the multiplets, including the S z - 
degeneracy, is 2 N . 

In the following we discuss the influence of the cor- 
relations between the electrons on the total transition 
probability 
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(4) 



where the spins S Z i and S z j refer to the states and 
(^jl, respectively. The matrix elements cj, a \^>i) im- 
ply first of all a spin selection rule, namely that each 
added or removed electron can change both the total 
spin S and the magnetic quantum number S z only by 
±1/2. In Jl(],[ll| these selection rules were included via 
the Clebsch-Gordan coefficients. They describe the com- 
bination of the initial spin state (Si, S Z i) with an electron 
(1/2, ±1/2) to form the final state (Sj,S Z j) and yield 



C 



3,1 



S i + 1 S{ 



25. 



Ss s ,Si-l/2 (5) 



after summation over S z i and S z j. This approach ig- 
nores the spatial degrees of freedom which make the Mji 
considerably different from the C'j.i- 

In tables || and III the numerically calculated Mjj and 
the corresponding Cj i are shown for transitions between 
two and three and three and four electrons in the dot, 
respectively. The truncation of the Hilbert space to M < 
13 leads to an absolute error of 10^3 (table H or 10 
(table III). In a few cases we checked the improvement 
achieved by using M=14 single electron levels. 

Some transitions are almost completely suppressed in 
the case of the Mj :i compared to the Cj t %. The two rea- 
sons can be seen from 



i*i> 



E b "* b ^ E^i^i"*) 



(6) 



Firstly it may be impossible to create (yA by adding 
one electron to \vi) (spatial selection rule) for the largest 
jb^by . |. Secondly, the various contributions to the sum- 
mation over Vi and Vj may cancel due to different signs 



2 



of h Vi h v . . This latter cancellation seems to have been 
neglected in ||. As an example, table IV shows the four 
main contributions for the 



l*§> 



/J 72 ) for N = 3 and 
for N — 4 together with their corresponding electronic 
configurations (see table | for the notation). They are 

Of the sixteen pos- 
only two (marked 



needed to calculate M^^ in table III. 
sible transitions shown in table 
with * and fulfill the spatial selection rule. In addition 
their contribution to M^^ cancel each other due to oppo- 
site signs (M 3)0 ~ (-0.35 x 0.29 + 0.33 x 0.39) 2 = 0.0006). 
This explains the smallness of the corresponding tran- 
sition probability. An other example is the transition 
between the first excited state for N — 3 and the ground- 
state for N = 2, (*f =1/2 |c t |*5 ; = ), table H In this case, 
we were able to follow the evolution of the many-electron 
states down to zero interaction, ending at Slater deter- 
minants \vi) and \vj) between which the spatial selection 
rule forbids transitions. 

The effect on non-linear transport properties is demon- 
strated in figure [l]. It shows grey scale plots of the dif- 
ferential conductance versus the gate voltage Vq and the 
transport voltage. In figure |l](left) the T^[ R were as- 



sumed to be proportional to C, 



while the calculated 



Mj i (eq. ||) were used in figure |l|(right). Grey areas cor- 
respond to regions of zero differential conductance. Black 
and white lines are related to positive and negative differ- 
ential conductances, respectively (spin blockade fl(||lT|| ). 
They reflect excited many electron states that become 
available for transport when gate and/or bias voltages 
are increased. On average, the number of lines is re- 
duced in figure |l|(right) as compared to figure 0(left) (eg. 
black arrows) . This reflects the suppression of transition 
matrix elements by the spatial selection rule discussed 
above. In some regions, however, the differential conduc- 
tance is even enhanced (cf. white arrows in figure ^). 
This is caused by considerable upheavals in the station- 
ary occupation probabilities Pj obtained from the rate 
equation when the full matrix elements are considered in 

(§■ 

The conductance peaks at low transport voltages 
shows different peak heights as presented in figure |^. 
This is directly related to the spatial properties of the 
many body states. Other works |^4] , |l5| explain this 
feature, also observed experimentally P,p[, within the 
framework of noninteracting electrons picture by semi- 
classical chaotic motions. In how far this picture can be 
generalized to the correlated electron situation deserves 
further research. Similar results in the presence of mag- 
netic fields were shown in |l6| . 

In summary, we have studied the electron transport 
through a quantum dot taking fully into account the cor- 
related eigenstates of the interacting electrons inside the 
dot. The spatial selection rule is shown to explain the 
suppression of certain transitions between N and (JV±1) 
electron states that would be allowed when taking into 



account only the spin selection rules. Despite the ob- 
tained tendency towards reduced transition probabilities 
Mj,i some of the peaks in the differential conductance 
are even enhanced as a result of considerable changes 
in the stationary occupation probabilities. Furthermore, 
the correlations between the electrons induced by the 
Coulomb interaction lead to characteristic variations in 
the heights of the linear conductance peaks. Our results 
show that non-linear transport spectroscopy provides in 
principle valuable information about the correlated dot 
states. To extract this information, and the correspond- 
ing physics, however, very careful theoretical modelling 
is required. 
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No.SCC*-CT90-0020 and CHRX-CT930126). 
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FIG. 1. Differential conductance versus transport and 
gate voltages (in units of e/ag) in linear grey scale (dark: 
positive; bright: negative). The electron number inside the 
diamond shaped Coulomb blockade region is N = 3. Left: 
transition probabilities proportional to C7j,». Right: transition 
probabilities proportional to Mj,j. The arrows are explained 
in the text. 
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FIG. 2. Current versus gate voltage for small transport 
voltage (V < A) using eq. (H). Temperature is WmK. The 
first peak corresponds to oscillations between N = 1 and 2 
electrons. 
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TABLE I. Sequence of increasing energy eigenvalues 7?f 
together with their total spins S. 
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TABLE II. Comparison between numerically calculated 
matrix elements Mj,i, eq.(Q), and corresponding values ob- 
tained by neglecting the spatial part of the wave function, 
Cj t i, eq. for the non-vanishing transition probabilities 

l*f(^=2)) - 1*^3))- 
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TABLE III. Same as table [o] for tran- 

sitions |*f(jv=3)) — * l*^^ 1 / 2 ))- Different columns are used 
for state with same spins S but different energies (see table 

0)- 
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TABLE IV. The four largest expansion coefficients 
^{jv,e ! ,s_} Q f |yy 2 ) anc l |$0) needed to calculate the entry 
0.01 in table |l|. The basis states \v) are illustrated for TV = 3, 
4 according to the occupations of single electron levels. Only 
two transitions between these states are possible by creating 
or annihilating one electron, marked with * or with . 
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